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Abstract 

We compute derived quantities for various values of the model parameter of the Cornell potential 
model for the S-wave heavy quarkonia with radial quantum numbers n = 1, 2, and 3. Our results 
can be used to determine leading and relative-order-?; 2 nonrelativistic quantum chromodynamics 
matrix elements for 5*-wave charmonia and bottomonia such as ip(2S), r/ c (2S), and T(nS) for 
n = 1, 2, and 3. These matrix elements will be essential ingredients for resumming relativistic 
corrections to processes involving those S- wave heavy quarkonium states. 

PACS numbers: 12.38.-t, 12.39.Pn, 14.40.Gx 
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I. INTRODUCTION 



One of the most interesting recent developments in heavy-quarkonium phenomenology 
is the introduction of a new technique for resumming relativistic corrections to S-wave 
quarkonium production and decay rates. The technique resums corrections to all orders in 
the heavy-quark velocity v in the heavy- quark- ant iquark (QQ) rest frame within the 

color-singlet mechanism of the nonrelativistic quantum chromodynamics (NRQCD) factor- 
ization approach [3j. With the method, the relative-order-f 2 NRQCD matrix element for the 
.S-wave charmonium, which has a power-ultraviolet divergence and needs subtraction, has 
been evaluated with improved accuracy [l|, compared with lattice calculations ^] and with 
the determination [5( obtained by using the Gremm-Kapustin relation 0. In Ref. fl, the 
generalized version of the Gremm-Kapustin relation within the Cornell potential model ?! 
was derived. The generalized Gremm-Kapustin relation allows one to resum a class of rela- 
tivistic corrections in a potential-model color-singlet QQ wave function. The resummation 
method has also been applied to determine leading-order NRQCD matrix elements for the 

tant values for the matrix elements are significantly greater 
2j. In addition, the method has provided a reasonable so- 
to the long-standing puzzle of the cross section for e + e~ — > J /ip + rj c measured 
at the .B-factories [9(, which has been one of the greatest discrepancies between theory and 
experiment within the Standard Model. 

Therefore, it is worthwhile to extend this method to study radially excited 5*-wave char- 
monia, ip(2S) and r] c (2S), and spin-triplet S-wave bottomonia, T(nS), with radial quantum 
numbers n = 1, 2, and 3. Unfortunately, the potential-model parameters and derived quan- 
tities reported in Ref. [l| are only for the IS* and the 25* states. The tabulation of the 
parameters in Ref. |l| is not convenient to use in combination with the improved version of 
the resummation method reported in Ref. [2], which considered only the IS charmonium 
states. 

In this paper, we compute derived quantities for various values of the model parameter 
of the Cornell potential model for S"-wave heavy quarkonia with radial quantum numbers 
n = 1, 2, and 3. The Cornell potential is a linear combination of the Coulomb and linear 
potentials. Following Ref. Q, we rescale the Schrodinger equation and solve the equation 
numerically for various values of the model parameter A, which determines the strength 



15* charmonium states. The resu 
than those known previously [l 
lution 
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of the Coulomb potential relative to the linear potential. To achieve improved accuracies 



in calculating energy eigenvalues, we use the numerical method 1 given in Refs. ll|, [12]. 
The potential-model parameters and derived quantities are listed as functions of A. These 
parameters are ready for use to determine leading and relative-order-f 2 color-singlet NRQCD 
matrix elements for both S'-wave charmonia and bottomonia. 

The remainder of this paper is organized as follows. In Sec. [TTJ we give a brief description 
of the Cornell potential model. Our numerical results are given in Sec. IHfl which is followed 
by a discussion in Sec. [IVl 



II. POTENTIAL MODEL 



In order to compute leading and higher-order color-singlet NRQCD matrix elements for 
an S'-wave heavy quarkonium with the radial quantum number n = 1, 2, or 3, we need to 
compute the binding energy e n s of the nS state that appears in the generalized Gremm- 



Kapustin relation 



ll, |6j . In this section, we describe briefly the Cornell potential model [?| 
that we use to compute e n g- We refer the reader to Refs. l|, |2[] which contain more complete 



descriptions of the model in conjunction with the NRQCD factorization formalism and the 
resummation technique for relativistic corrections to quarkonium processes. 

We employ the Cornell potential [7], which parametrizes the QQ potential V(r) as a 
linear combination of the Coulomb and linear potentials: 



V(r) 



K 

r 



ar, 



(1) 



where k is a model parameter for the Coulomb strength and a is the string tension. The 
relation between the string tension a and the corresponding parameter a in the original 
formulation of the Cornell potential model Q is a = \j\fo. By varying the parameters in 
the Cornell potential, one can obtain good fits to lattice measurements of the QQ static 
potential 13J]. Therefore, we assume that the use of the Cornell parametrization of the QQ 
potential should result in errors that are much less than the order-w 2 errors (about 30% for 
a charmonium and about 10% for a bottomonium) that are inherent in the leading-potential 
approximation to NRQCD [l, 2|. 



1 In previous calculations in Refs. [H, the authors quoted the energy eigenvalues determined in Ref. 
A recently improved analysis given in Ref. |2] employed the numerical method used in the present work. 



The Schrodinger equation for the S"-wave radial wave function R n s(r) with the radial 
quantum number n is 



mr 2 dr \ dr 



R n s(r) = e nS R n s(r), 



(2) 



where m is the quark mass and e n s is the binding energy for the nS state. The model param- 
eter m is distinguished from the heavy-quark mass tuq that appears in the short- distance 
coefficients of NRQCD factorization formulas. For an S"-wave state, the wave function is 
^ nS (r) = Rns(r)/V^r. 

The Schrodinger equation in Eq. ([2]) depends on the model parameters m and k, where we 
assume that the string tension a is common to both charmonium and bottomonium states 
and that the value for a can be determined from lattice measurements of the QQ static 
potential. Then, the dependence on the flavor appears through m and n. Introducing the 
scaled radius p and scaled coupling A 3], 



p = (<rm) 1//3 r, 
A = " 



( a / m 2)l/3' 

which are dimensionless, one can rewrite the radial equation in Eq. ([2]) as [3] 

d 2 A 



(3a) 
(3b) 



-j-2 H P + CnS 

dp z p 



Uns(p) = 0, 



(4) 



where u n s{p) and ( n s are the dimensionless radial wave function and the dimensionless 
energy eigenvalue of the nS state, respectively. The relation between R n s( r ) and u n s(p) is 

Uns{p) 



Rnsir) 



'am ■ 



P 



(5) 



where the wave functions are normalized according to 



Ks(p)\ 2 dp= I \R n *(r)\ 2 r 2 dr = 1. 



'0 JO 

The binding energy is related to the dimensionless eigenvalue ( n $ as 



(6) 



s nS = [a 2 /m] l /\ nS {\). 



(7) 



Note that Eq. (j4j) depends only on A. Therefore, the scaled equation can be solved for 
a given A to get the wave function u n s(p) and the eigenvalue ( n s- The flavor dependence 
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appears when we invert them to get the radial wave function R n s(r) and the energy eigen- 
value e n s- In this step, m looks independent of the model parameter A. However, we can 
express m in terms of a, A, and the 1S-2S mass splitting [l|, |2j, |7|: 

"C 2 s(A) - Gs(A) n:! 



m(A) 



<j 



(8) 



m 2S - mis 

For S'-wave states, the wave function at the origin ip n s(fy = R n s(0) / \J can be expressed 
as [ll,y,LZ| 



m 



hM0)| 2 = - / d^ nS (r)\ 



,dV(r) crm(X) 



dr 



Air 



[1 + XF nS (X)} 



(9) 



where F n s(X) is the expectation value of 1/p 2 for the nS state: 

dp 



Fns(X) 



P 



Ks(p)T 



(10) 



For purposes of computation of the NRQCD matrix elements, it is convenient to express 
those matrix elements in terms of the potential-model parameters listed above. A convenient 
parametrization can be found in Ref. [2]. According to Ref. |2|, the leading-order NRQCD 
matrix element for the S'-wave heavy quarkonium depends on A, m(X), and F n s(X). The 
ratio of the relative-order-f 2 NRQCD matrix element to the leading-order one is propor- 
tional to Cns(A). The electromagnetic decay rate for the quarkonium, in which relativistic 
corrections to all orders in v are resummed, is expressed in terms of these NRQCD matrix 
elements [2]. Because m(A), F nS (X), and ( n s{X) depend on A, the decay rate is completely 
determined by A. Determination of the best value for A can be made by imposing a re- 
quirement that the resummed formula for the decay rate should reproduce the measured 
rate. Therefore, the leading-order and the relative-order-f 2 NRQCD matrix elements are 
completely determined [2]. 



III. NUMERICAL RESULTS 



In this section, we list our numerical values for the parameters and derived quantities 
for the S- wave heavy quarkonia with radial quantum numbers n = 1, 2, and 3. In order to 
provide a set of parameters that can be used for both charmonia and bottomonia, we list the 
values for ( n s and F n $ that are common to both cases. In Table [H we tabulate those values 
as functions of the model parameter A, whose range has been chosen so that the parameters 

5 



TABLE I: Scaled energy eigenvalues CnS and F n g of the 5-wave heavy quarkonium with radial 
quantum numbers n = 1, 2, and 3 as functions of the Coulomb strength parameter A. 



A 


Cis 


Fis 


(2s 


F 2S 


C3S 


F 3S 


0.0 


2.338107 


1.1248 


4.087949 


0.8237 


5.520560 


0.6983 


0.1 


2.253678 


1.1869 


4.029425 


0.8525 


5.473169 


0.7178 


0.2 


2.167316 


1.2532 


3.970286 


0.8821 


5.425462 


0.7375 


0.3 


2.078949 


1.3237 


3.910531 


0.9125 


5.377441 


0.7576 


0.4 


1.988504 


1.3989 


3.850160 


0.9437 


5.329112 


0.7780 


0.5 


1.895904 


1.4789 


3.789174 


0.9756 


5.280478 


0.7987 


0.6 


1.801074 


1.5641 


3.727575 


1.0083 


5.231545 


0.8196 


0.7 


1.703935 


1.6546 


3.665364 


1.0417 


5.182316 


0.8408 


0.8 


1.604409 


1.7507 


3.602543 


1.0758 


5.132798 


0.8622 


0.9 


1.502415 


1.8527 


3.539116 


1.1105 


5.082996 


0.8837 


1.0 


1.397876 


1.9608 


3.475087 


1.1459 


5.032914 


0.9054 


1.1 


1.290709 


2.0753 


3.410458 


1.1818 


4.982560 


0.9272 


1.2 


1.180834 


2.1965 


3.345233 


1.2183 


4.931938 


0.9492 


1.3 


1.068171 


2.3246 


3.279418 


1.2552 


4.881053 


0.9712 


1.4 


0.952640 


2.4599 


3.213016 


1.2927 


4.829913 


0.9933 


1.5 


0.834162 


2.6026 


3.146031 


1.3306 


4.778522 


1.0155 


1.6 


0.712658 


2.7529 


3.078468 


1.3689 


4.726886 


1.0377 


1.7 


0.588049 


2.9110 


3.010330 


1.4077 


4.675010 


1.0598 


1.8 


0.460260 


3.0773 


2.941621 


1.4468 


4.622899 


1.0820 


1.9 


0.329215 


3.2518 


2.872344 


1.4862 


4.570560 


1.1041 


2.0 


0.194841 


3.4348 


2.802503 


1.5260 


4.517996 


1.1262 


2.1 


0.057065 


3.6264 


2.732099 


1.5662 


4.465212 


1.1482 


2.2 


-0.084182 


3.8269 


2.661134 


1.6066 


4.412212 


1.1702 


2.3 


-0.228969 


4.0364 


2.589611 


1.6474 


4.359001 


1.1921 


2.4 


-0.377362 


4.2550 


2.517529 


1.6885 


4.305582 


1.2139 


2.5 


-0.529425 


4.4829 


2.444888 


1.7299 


4.251959 


1.2356 


2.6 


-0.685221 


4.7202 


2.371688 


1.7716 


4.198135 


1.2573 


2.7 


-0.844808 


4.9670 


2.297928 


1.8137 


4.144112 


1.2788 


2.8 


-1.008244 


5.2235 


2.223605 


1.8561 


4.089893 


1.3002 


2.9 


-1.175584 


5.4897 


2.148718 


1.8989 


4.035481 


1.3216 


3.0 


-1.346882 


5.7657 


2.073261 


1.9421 


3.980877 


1.3428 
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can be used to determine the NRQCD matrix elements for both charmonia and bottomonia. 
The dependence of the eigenvalues CnS and F n s on A are shown in Figs. [Hand [21 respectively. 
In Ref. the authors solved the Schrodinger equation in Eq. numerically by using 



the inverse iteration method introduced in Ref. [12] . The method uses a trial wave packet 
which is a linear combination of the eigenfunctions for a given system and the amplification 
operator of the form (H — Ctriai) -1 , where H is the Hamiltonian of the system and Atrial 
is a trial eigenvalue. If this operator acts on an eigenfunction of the Hamiltonian, then 
the operator yields a multiplicative factor (£ — Ctriai) 1 , where C is the eigenvalue of that 
eigenfunction. As Atrial gets closer to C, the multiplicative factor blows up. Therefore, if 
this operator acts on a wave packet, then an eigenfunction, whose eigenvalue is the nearest 
to Ctriai, is selectively amplified. Therefore, by applying the operator repeatedly, one could 
obtain the eigenfunction and the corresponding eigenvalue from the expectation value of the 
Hamiltonian. By varying Ctriai, one could also obtain the energy spectrum. The details of 
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121 ] . Here, we use this method to generate 



the numerical method can be found in Refs. 
TableE 

If one substitutes ( n s{ty listed in Table Hand the measured 1S-2S mass splitting for either 
the charmonium or the bottomonium into Eq. (jSJ), one can obtain the mass parameter m(A) 
as a function of A. e n s and ip n s(fy are calculable by substituting the parameters Cns(A), 



m(A), and F nS (X) into Eqs. (J7J) and Qj, respective 
optimal value for the parameter A is given in Refs. 



v. A convenient way to determine the 

m 



IV. DISCUSSION 



We have computed derived quantities for various values of the model parameter of the 
Cornell potential model for the ,5-wave heavy quarkonia with radial quantum numbers n — 1, 
2, and 3. The scaled Schrodinger equation is solved numerically for various values for the 
model parameter A which determines the strength of the Coulomb potential relative to the 
linear potential. The scaled energy eigenvalue CnS and a derived potential quantity F nS 
are computed as functions of A. These numbers are useful in determining the leading and 
relative-order-f 2 NRQCD matrix elements for both ^-wave charmonia and bottomonia with 
radial quantum numbers n = 1, 2, and 3. As an application, the leading-order NRQCD 
matrix elements for the T(nS) have already been calculated and used to determine the 
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FIG. 1: Eigenvalues (, n s as functions of A. 



branching fractions and charm-momentum distributions for the inclusive charm production 
in T(nS) decays, which are being analyzed by the CLEO Collaboration 141 ] . The result 
listed in this paper can also be used to determine the NRQCD matrix elements for the 2S 



charmonium states. Once recent studies 



15 



16 



171 ] to compute relativistic corrections 



to the leptonic width of the S'-wave spin-triplet quarkonium are extended to complete the 
resummation of relativistic corrections to the process at next-to- leading order in a s , our 
results can also be used to determine the leading-order and relative-order-f 2 NRQCD matrix 
elements for the S'-wave quarkonium with accuracies better than the best available values 
in Ref. Q. 
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FIG. 2: F n s as functions of A. 
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